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Abstract 

Using the full data set of the NA48/2 experiment, the decay K ± — > ■K ± e + e^j is 
observed for the first time, selecting 120 candidates with 7.3±1.7 estimated background 
events. With — * ^ir^ as normalisation channel, the branching ratio is determined 
in a model-independent way to be Br(X ± — > ir ± e + e~j, m ee7 > 260MeV/c 2 ) = (1.19± 
0.12 stat ± 0.04 syst ) x 10~ 8 . This measured value and the spectrum of the e + e~7 
invariant mass allow a comparison with predictions of Chiral Perturbation Theory. 
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1 Introduction 



The decay K ± — > 7r ± e + e 7 is similar to the decay — ► 7r 77, with one of the pho- 
tons internally converting into a pair of electrons. Both decays can be described in the 
framework of Chiral Perturbation Theory (ChPT). The lowest order terms are of order 
p A , where predominantly loop diagrams contribute to the amplitude [1]. This leads to a 
characteristic signature in the e + e~7 invariant mass distribution, which is favored to be 
above 2m n + and exhibits a cusp at the 2m n + threshold. In ChPT, the loop contribution 
is fixed up to a free parameter c, which is a function of several strong and weak coupling 
constants and expected to be of 0(1) [2]. 

Higher order ChPT calculations on K + — > 7r + 77 have been performed, but are model- 
dependent. Theoretical predictions for K + — ► 7r + e + e~7 exist [3], following the earlier 
work on K + — ► 7r + 77 [2]. The predicted branching ratios lie in the range between 0.9 and 
1.7 x 10 -8 , for values of |c| < 2. Experimental results are available only for K + — > 7r + 77, 
based on the observation of 31 signal candidates by the E787 experiment [4]. 

In this letter, we report the first observation of the decay — > 7r ± e + e _ 7 and the model- 
independent measurement of its branching fraction, using 

K ± ^ n ± n D 

with 7T£, — > e + e 7 

as the normalisation channel. These results have been derived from the full data set of the 
NA48/2 experiment. 

2 Experimental Set-up 

The NA48/2 experiment took data in 2003 and 2004 at the CERN SPS. Two beams of 
charged particles were produced by a 400 GeV/c proton beam impinging on a Be target in 
a 4.8 s long pulse repeated every 16.8 s. Positive and negative particles with momenta of 
(60 ±3) GeV/c were simultaneously selected by an achromatic system, which split the two 
beams in the vertical plane and then recombined them on a common axis. After passing 
through a collimator, the beams were split and recombined again in a second achromat. 
Finally, the two beams passed a cleaning and a defining collimator before entering the decay 
volume housed in a 114 m long evacuated tank with a diameter between 1.92 and 2.4 m 
and terminated by a 0.3% radiation lengths thick Kevlar window. The axes of both beams 
coincided within 1 mm inside the decay volume. The beams were primarily composed of 
charged pions, with a fraction of 5 — 6% of . On average, about 4.8 x 10 5 K + and 
2.7 x 10 5 K~ per pulse decayed in the fiducial decay volume. A more detailed description 
of the beamline can be found in [5]. 

The decay region was followed by the NA48 detector [6]. The momenta and positions 
of charged particles were measured in a magnetic spectrometer. The spectrometer was 
housed in a helium gas volume and consisted of two pairs of drift chambers before and 
after a dipole magnet with vertical magnetic field direction, giving a horizontal transverse 
momentum kick of 120 MeV/c. Each chamber had four views (x, y, u, v) with two sense 
wire planes in each view. The u and v views were inclined by ±45° with respect to the x-y 
plane. The space points, reconstructed by each chamber, had a resolution of 150 /im in each 
projection. The momentum resolution of the spectrometer in 2003/2004 was measured to 
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be (Jp/p = 1.02% © 0.044% x p } with p in GeV/c. The magnetic spectrometer was followed 
by a segmented plastic scintillator hodoscope with one plane of vertical and one plane of 
horizontal strips, respectively. It was used to produce fast trigger signals and to provide 
precise time measurements of charged particles. The time resolution of the hodoscope was 
better than 200 ps. 

Photon and electron energies were measured with a 27 radiation length thick liquid- 
krypton electromagnetic calorimeter (LKr). It was read out longitudinally in 13248 cells of 
2x2 cm 2 cross-section. The energy resolution was determined to be cte/E = 3.2%/ 'y/E © 
9%/i?ffi0.42%, with E in GeV. The spatial and time resolutions were better than 1.3 mm 
and 300 ps, respectively, for photon and electron clusters above 20 GeV. 

Additional detector elements, such as the hadron calorimeter and the muon and photon 
veto counters, were not used in the present analysis. 

The events were selected by a two-level trigger which was optimised for events with 
three charged tracks. At the first level, three-track events were triggered by requiring 
coincidences of hits in the two hodoscope planes. The second level trigger was based on 
the hit coordinates in the drift chambers. It required at least two tracks to originate from 
the decay volume with a reconstructed distance of closest approach of less than 5 cm. 
The trigger efficiency for the selection of the normalisation channel — > tt^tt^ was 
(96.48 ± 0.05)%, determined from data events taken with a down-scaled complementary 
trigger. 

In total, NA48/2 collected about 18 x 10 9 triggers. In the course of data-taking, the 
magnet polarities of both the beamline and the spectrometer were regularly reversed to 
have similar conditions for decays of positive and negative kaons. 

3 Monte Carlo Simulation 

In order to compute the acceptance of signal, normalisation, and background channels, a 
detailed GEANT-based [7] Monte Carlo simulation was employed, which included the full 
detector and material description, stray magnetic fields, drift chamber inefficiencies and 
misalignment, and beamline simulation. 

For the signal channel, the full matrix element was used [3], with a value of c = 1.8, in 
agreement with the measurement of K + — > 7r + 77 [4]. 

For all other channels, if not otherwise mentioned below, the known theoretical ma- 
trix elements were used. Radiative corrections were applied to the simulation of signal, 
normalisation, and A'- — > 7r ± e + p _ by using the PHOTOS package [8j. 

4 Data Selection 

The analysis described here is based on the full data set of the NA48/2 experiment, recorded 
in 2003 and 2004. The selection of the signal events was performed in two steps. At first, 
described in the next section, a set of basic selection criteria was applied to define the signal 
region and to assure the quality of the selected candidate events. In a second step, described 
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in the following section, further selection criteria were applied to suppress contributions of 
the various background sources. 

4.1 Event Selection 

Each selected event had to have at least one combination of three tracks with a total charge 
of ±1 and one cluster in the LKr calorimeter not associated with any track. Each track 
was required to have a radial distance from the detector axis of at least 12 cm in the first 
drift chamber and to lie well inside the active LKr calorimeter region, i.e. well inside its 
outer edge and more than 15 cm from the detector axis. The distance between any two 
tracks in the first drift chamber had to exceed 2 cm. This latter requirement rejects all 
events with external photon conversions in the detector material before the spectrometer. 

The three tracks had to be compatible with originating from the same decay vertex and 
to have a distance of closest approach of less than 4 cm for each of the three pairs of two 
tracks. The longitudinal position of the reconstructed decay vertex had to be more than 
2 m and less than 98 m down-stream of the final collimator and within a radius of 3 cm 
around the beam axis. The track times, measured in the scintillator hodoscope, had to 
be at most ±3 ns from the mean of the track times. In less than 1% of the events, no 
hodoscope information was available for at least one track. For those events, the track 
times measured in the drift chambers were taken and required to be at most ±6 ns from 
the mean track time. 

Pions and electrons were identified by the ratio E/p of energy deposited in the LKr 
calorimeter and momentum measured in the spectrometer. Electrons and pions were re- 
quired to have E/p > 0.94 and E/p < 0.8, respectively. With these requirements the 
probability for mis-identification of electrons or pions is of the order of a few per mille. 
Two tracks of opposite charge had to be identified as electrons with each having a momen- 
tum greater than 3 GeV/c. The third track had to be a pion candidate and had to have a 
momentum greater than 4 GeV/c. 

Photon candidates were defined as calorimeter clusters unassociated to charged tracks 
and required to lie 15 cm from the detector axis and well inside the outer edge of the LKr. 
The distance to the projected impact point of any pion candidate had to exceed 25 cm, 
the distance to the electron tracks or any other possible cluster had to exceed 10 cm. The 
reconstructed cluster energy had to be greater than 3 GeV, and the time difference to the 
mean of the cluster time and the track times measured in the drift chambers had to be 
smaller than 6 ns. 

The sum of pion, electron, positron, and photon momenta was required to be between 54 
and 66 GeV/c. An energy centre-of-gravity (x cog , y cog ) = %iEi / J2i E i> Y,i Vi E i I J2i E i) 
was defined by using the transverse positions xt and yi and the energies Ei of the projected 
tracks and the photon cluster at the front surface of the LKr calorimeter. The tracks were 
projected onto the LKr surface from their positions and directions in the first drift chamber 
before the spectrometer magnet. The radial distance between the energy centre-of-gravity 
and the beam had to be less than 3 cm. 

To suppress background events coming from decays with more than one photon, we re- 
quired that no other unassociated cluster was in-time with the event. Since this would reject 
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also events with photons from bremsstrahhmg on detector material or with shower fluctua- 
tions of pion showers, we still allowed additional clusters with E [GeV] < 7 — 0.14d e [cm] or 
E [GeV] < 15 — 0.25 [cm], where d e and d v are the distances of the cluster to the impact 
point of an electron and the pion track, respectively. This requirement against additional 
unassociated clusters rejected about 0.3% of all events. 

With these basic selection criteria, a sample of about 22.8 million events was obtained. 
At this level of the selection, the data were dominated by — ► tt^tt^ decays, where the 
7T° underwent a Dalitz decay -k q d — ► e + e~j. 

4.2 Background Suppression 

A number of additional selection criteria had to be applied to effectively suppress the 
remaining background. In the following, we examine the possible sources of background 
to the signal. 

• K± — 71^71-°,, K ± -» ir^e^v, and K± -> n^fj^u decays 

The decay — ► tt^tt^ {K^d) with tt° d — > e + e~j has exactly the same signature as 
the signal channel. We therefore rejected events for which 120 MeV/c 2 < m e + e - 7 < 
150 MeV/c 2 . To evaluate this requirement, we assigned the electron mass to each 
track and applied the cut to both opposite-charged track combinations. This com- 
pletely rejects also the semileptonic decays — ► ix^e^v and A'± —> TT^fi^u as well 
as the small amount of doubly-misidentified K^d events, where both the pion and 
an electron are misidentifled. 

• 77^7707 decays 

The decay — > 7r ± 7r£,7 consists of two amplitudes: Inner Bremsstrahlung (IB) 
and Direct Emission (DE). If the radiative photon is lost, the decay is rejected by 
the cut against — > ^7r° D decays. However, if the photon of the tt q d decay is lost, 
the decay may fake a signal event. — > tt^tt^ events are the major background 
source and contribute with 3.1 ± 0.5 IB and 0.12 ± 0.03 DE events, as determined 
from the simulation, using the measured rates of IB and DE transitions [9]. 

• K ' — ¥ 7r = ' = 7r e+e — decays 

The decay J fC ± — > 7r =l: 7r e + e , which comes from K — > 7r ± 7r°7 with an internal 
conversion of the additional photon, has not been measured yet. By evaluating the 
internal conversion probability, we estimated its branching fraction to be half of that 
of — > 7r ± 7T£ ) 7. Due to the uncertainty of the estimation, we assigned a ±50% 
systematic uncertainty to this value. — > ir ± Tr°e + e~ events were simulated by 
modifying the — > 7r ± 7r°7 simulation; the conversion of the photon was added by 
generating the photon mass with a probability density proportional to the inverse 
square of the photon mass. The ratio of IB and DE amplitudes was taken from 
— ► 7r ± 7r°7 decays [9]. From this, we estimated the amount of background from 
— > 7r ± 7r°e + e to 1.6 ±0.5 ±0.7 events, where the second uncertainty comes from 
the estimation of the branching fraction. The contribution from DE is practically 
negligible. 
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Smaller angle(e,y) - Angle(e,e) [rad] 



Figure 1: Difference of the smaller of the angles between the photon and e and the angle 
between e + and e~ in the e + e~7 rest frame for data (circles) and signal and 
radiative — > ir ± e + e~ MC simulation. 

• Radiative K^ 1 —* 7v^e + e~ decays 

The probability of the rare decay K — > ir ± e + e~ with an observable 7 from internal 
or external bremsstrahlung is of 0(1CT 8 ), similar to the signal channel. To reject 
these events, we made use of the different decay kinematics in the e + e~7 rest frame. 
For the signal, the photon repels from the e + e~ system, while in case of — ► 
7r ± e + e~ + 7brems the electrons, in the e + e~7 rest frame, fly back-to-back, and the 
photon is close to one of them. We therefore required <£(e + ,e~) < min(<j(e , 7)) 
for the angles between e + and e~ and between e ± and the photon, respectively, 
in the e + e~7 rest frame (see Figure [1]). The remaining background from — > 
7r ± e + e~ decays was determined to be 0.8 ± 0.5 events from MC simulation. This 
estimate includes a systematic error of ±50%, which reflects a disagreement between 
— y 7r ± e + e~ data and MC in event numbers in the region with <(e + ,e~) > 
min(<(e ± , 7)). 

• K± -> t^tt ^ decays 

The decay — ► 7r ± 7r°7r^ ) was strongly suppressed by the cut on additional clusters 
and the rejection of the it° D decays. Its contribution to the signal region was estimated 
by MC simulation to be 0.7 ± 0.7 events. 

• Accidental activity 

Accidental overlap of separate events may fake signal events. To estimate the amount 
of such events in the signal sample, we studied the sidebands of the time distributions 
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in the hodoscope and the LKr calorimeter. One event was found in the calorimeter 
time sideband, which corresponds to a background estimation of 1 ± 1 events from 
accidental activity. 

Other potential sources of background as e.g. — > 7r ± 7r + 7r~(7), — > 7r + 7r~ e^v(y), 
or — > 7r°7r^ ) e ± i/(^) were found to be irrelevant. 

The signal region was defined by requiring 480 < m 7r ± e + e - 7 < 505 MeV/c 2 . Since ChPT 
predicts only small signal rate and the background increases for low values of m e + e - 7 , 
we also required "i e + e - 7 > 260 MeV/c 2 . 120 signal candidates were found, including an 
estimated total background of 7.3 ±1.7 events. The background channels are listed in 
Table Q] together with their respective branching fractions and expected contributions to 
the signal. All background expectations were obtained by normalising to the total kaon 
flux. 

The signal acceptance, as determined from MC simulation, depends on m e + e - 7 . It was 
between 6% and 7% for 260 < m e + e - 7 < 330 MeV/c 2 , and decreased to 2.5% for events 
near the kinematical edge. The projections of the signal candidates on m 7r ± e + e - 7 and 
m e + e - 7 , together with the background contributions, are shown in Figure [2j 



Background source Branching ratio Expected events 



K± -> 7 r ± 7r^7 (IB) 


3.3 x lO" 6 


3.1 ±0.5 


7r ± 7r2,7 (DE) 


5.3 x 10~ 8 


0.12 ±0.03 




~ 1.7 x 10" 6 


1.6 ±0.9 


-» 7r ± 7r°e + e- (DE) 


~ 2.6 x 10~ 8 


0.02 ±0.01 




2.9 x 10~ 7 


0.8 ±0.5 


K ± ^ ^o^o 


2.1 x 10~ 4 


0.7 ±0.7 


Accidentals 




1.0 ±1.0 



Sum 7.3 ± 1.7 



Table 1: Relevant background sources and number of expected events in the signal region 
with all selection criteria applied. Except for — > 7r ± 7r e + e _ , the branching 
fractions are taken from [9]; for — ► -K^it ^ and — > 7r ± 7r°e + e _ they are 
defined for T w = 55 - 90 MeV. 

4.3 Normalisation Channel 

For the normalisation channel — > tt^tt^ exactly the same selection as for the signal 
was applied, but without the criteria on the e + e~7 invariant mass and the e + e~7 decay 
angles. Instead, the e + e~7 invariant mass was required to be within m^o — 35 MeV/c 2 
and m n o + 30 MeV/c 2 . The asymmetry of this cut takes into account the radiative tail 
in the m e + e - 7 distribution. We found about 18.7 million K2 n D candidates including an 
estimated background of about 0.5%. The acceptance of the selection was determined to 
be 5.0% from MC simulation. The branching fraction of the normalisation is Br(A' ± — > 
^71^,71% -> e+e~7) = (2.51 ± 0.07) x 10~ 3 [9j. From this, we determined the total flux 
of kaon decays in the fiducial volume to be &k = (1-48 ± 0.04) x 10 11 . 
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• Data 
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_ K + -> Ji + Ji0y 

□ K + ^> Ji + ji°e + e- 
I | K + ^> 7i + e + e _ 




0.45 0.5 0.55 

Invariant ^e+e y mass [GeV/c 2 ] 



^20 
o 
> 

0) 



| K + ^> 7t + 7t£>Y 

7i + 7t°e + e" 
7t + e + e" 



0.2 




0.25 0.3 0.35 

Invariant e + ey mass [GeV/c 2 ] 



Figure 2: Selected signal candidates and background expectation from MC simulation: 
(top) 7r ± e + e~7 invariant mass, (bottom) e + e~j invariant mass. The vertical 
lines indicate the accepted region for the branching ratio measurement. 



5 Results 



To determine the branching fraction in a model-independent way, we computed a partial 
branching fraction for each 5 MeV/c 2 wide m e + e - 1 interval i from 

jV 7ree7 - N bkg 1 



r7 ree7 n , » r bk, 



with the numbers N i and g of observed signal and estimated backgrounds events, 
and the signal acceptance A* eei in bin i. The overall trigger efficiency is e and <3?x the 
total kaon flux. By summing over the bins above "i e + e - 7 = 260 MeV/c 2 , we obtained 
B^K* -> 7r ± e + e-7,m e+e - 7 > 260 MeV/c 2 ) = (1.19 ± 0.12 stat ) x 10~ 8 , where the error is 
from data statistics only. This result is independent of the value of c and any theoretical 
assumption of the m e + e - 7 distribution. 

Several potential sources of systematic errors can affect the result and have been studied. 
The background estimation has a total uncertainty of ±1.7 events, as explained before, 
which results in an uncertainty of ±1.5% on the result. 

Possible imperfections of the description of the detector acceptance in the Monte Carlo 
simulation might also cause systematic effects on the branching fraction measurement. For 
an estimation of such effects we have varied the main selection cuts. To not fall victim of 
statistical fluctuations in the signal channel, the variations have only been performed in the 
normalization channel. This leads to a conservative estimate, since detector systematics 
are expected to cancel between signal and normalization. We found maximum changes of 
the result of the order of ±0.4%, which we assigned as the systematic uncertainty due to 
the detector acceptance. 

The particle identification via the E/p ratio could not be perfectly modelled in the 
simulation. However, the inefficiencies are expected to almost completely cancel between 
signal and normalization mode. The residual uncertainty on the electron identification was 
measured from K® 3 decays to better than 0.1%. The uncertainty of the pion identification 
efficiency was determined from variations of the E/p criterion in the normalization channel 
to be at most ±0.3%. Combining both, we assigned a total uncertainty of ±0.4% due to 
particle identification. 

The overall trigger efficiency should be the same for signal and normalisation to a great 
extent. A difference could only arise from the slightly different event topologies. Due to 
the lack of statistics, the trigger efficiency could not be measured for signal events. We 
therefore studied the dependency of the trigger efficiency of — > ^7^% events as a 
function of the event topology. From this, we obtained a systematic uncertainty of ±0.6%. 

The statistical error of the signal and normalisation MC samples contributes to ±0.9%. 

Finally, the external inputs of Br^K^ — > tt^tt^) and Br(7r£> — ► e + e~j) add an uncer- 
tainty of ±2.7% [9]. This is identical with the error quoted on the kaon flux in Section [43l 

All uncertainties of the measurement are listed in Table [2l The final result on the 
branching ratio is 

Br (IT* -» 7r ± e + e - 7, m e+e - 7 > 260 MeV/c 2 ) = (1.19 ± 0.12 stat ± 0.04 syst ) x 10~ 8 . 
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Source 


ABr/Br 


ABr [10" s ] 


Background subtraction 


± 1.5% 


±0.017 


Electron/pion identification 


± 0.4% 


± 0.005 


Detector acceptance 


± 0.4% 


± 0.005 


Trigger efficiency 


± 0.6% 


± 0.007 


MC statistics 


± 0.9% 


±0.011 


Normalisation 


± 2.7% 


± 0.032 


Total systematic uncertainty 


± 3.3% 


±0.04 


Statistical uncertainty 


± 9.7% 


±0.12 



Table 2: Summary of uncertainties of the branching ratio measurement. 

The distribution of the partial branching fractions is shown in Figure [3] and tabulated in 
Table 02 The quoted errors are confidence intervals for the unknown true value. We chose 
Pearson's \ 2 intervals for Poisson statistics [10] for them, defined as a± = \/ni + l/4± 1/2 
for each data bin with entry 7ij, before background subtraction, acceptance correction, and 
normalisation. The uncertainties on the background and acceptance estimates were added 
in quadrature in each bin, while the global systematic uncertainties — dominated by the 
normalisation — are not quoted in Figure [3] and Table 

We used the measured branching fraction and the shape of the e + e~7 spectrum to extract 
a value for the parameter c. Performing a least squares fit of the absolute prediction given 
in ref. [3] to the data with m e + e - 7 > 260 MeV/c 2 , we obtained c = 0.90 ± 0.45, where the 
error is dominated by the data statistics. The quality of the fit was x 2 / n doi — 8.1/17. This 
result is in agreement within about 1.2 standard deviations with the value of 1.8 ± 0.6, 
previously measured in K + — > 7r + 77 |4], and has a somewhat smaller error. Figure [3] shows 
the predicted spectrum for our best fit value of c and the previously found value, together 
with the background and acceptance corrected data. 

Using our measured value of c and ref. [3], we computed the differential branching fraction 
for m e + e - 7 < 260 MeV/c 2 and added it to our measured result. We then obtained for the 
total branching fraction Br(^ ± -> vi^e+e^) = (1.29±0.13 exp ±0.03 6 ) x 10 -8 , where the 
first uncertainty is the combined statistical and systematic error, and the second reflects 
the uncertainty in c. 
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Figure 3: Partial if ± — ► 7r ± e + e 7 branching fractions as function of the e + e 7 invariant 



mass. The signal region is defined for m e + e - 1 > 260 MeV/c 2 . The lines are 
the expectations from ref. [3] for c = 0.90 (solid, best fit) and c = 1.8 (dashed, 
estimate from K + — > 7r + 77 [4]). 





+ e - 7 interval 


A Br 


260 


- 265 


MeV/c 2 


0.17 


265 


- 270 


MeV/c 2 


0.38 


270 


- 275 


MeV/c 2 


0.40 


275 


- 285 


MeV/c 2 


0.63 


280 


- 280 


MeV/c 2 


0.62 


285 


- 295 


MeV/c 2 


0.94 


290 


- 290 


MeV/c 2 


0.95 


295 


- 300 


MeV/c 2 


0.95 


300 


- 305 


MeV/c 2 


1.07 



m e + e -^ interval A Br [10 9 ] 



+0.23 
-0.13 
+0.28 
-0.18 
+0.27 
-0.18 
+0.31 
-0.21 

+0.30 
-0.21 

+0.36 
-0.26 
+0.36 
-0.26 
+0.36 
-0.26 
+0.38 
-0.28 



305 
310 
315 
320 
325 
330 
335 
340 
345 



310 MeV/c 2 
315 MeV/c 2 
320 MeV/c 2 
325 MeV/c 2 
330 MeV/c 2 
335 MeV/c 2 
340 MeV/c 2 
345 MeV/c 2 
350 MeV/c 2 



1.09 
0.77 
0.68 
0.67 
0.32 
0.45 
0.89 
0.62 
0.25 



-0.39 

-0.29 
-0.34 
-0.24 
-0.33 
-0.23 
-0.35 
-0.25 
-0.25 
-0.14 

-0.29 
-0.18 
-0.41 
-0.28 
h0.40 
-0.25 
-0.43 
-0.17 



Table 3: Partial K± — > 7T e + e 7 branching fractions dependent on the e + e 7 invariant 
mass. Quoted are uncertainties from data and MC statistics and background esti- 
mation. All other uncertainties are completely correlated and amount to ±3.2%, 
dominated by the normalisation. 
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